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We show how the relativistic matter and velocity power spectra behave in different gauges. We 
construct a new gauge where both spectra coincide with Newtonian theory on all scales. However, in 
this gauge there are geometric quantities present which do not exist in Newtonian theory, for example 
the local variation of the Hubble parameter. Comparing this quantity to second order Newtonian 
quantities, we find that Newtonian theory is inaccurate on scales larger than 10 Mpc. This stresses 
the importance of relativistic corrections to Newtonian cosmological N-body simulations on these 
scales. 
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I. INTRODUCTION 

At early times the Universe was very close to isotropic 
and homogeneous, the best indication being the tiny 
fluctuations (about ~ 10"^) in the cosmic microwave 
background (CMB) radiation [l|. However, the struc- 
ture of the Universe that we observe today is very in- 
homogeneous on scales below the homogeneity scale of 
- 100 Mpc [113. There are galaxies (with mass collec- 
tion radius up to ~ 1 Mpc) and clusters of galaxies (at 
scales of ~ 10 Mpc). Up to scales of ^ 100 Mpc we find 
superclusters and large voids 

The theory of structure formation connects the early 
homogeneous Universe with the inhomogeneous Universe 
wc observe today. It is assumed that there exist initially 
small density perturbations, which grow due to gravi- 
tational attraction. The mathematical tool is relativis- 
tic cosmological perturbation theory [5|, which consid- 
ers small perturbations on a homogeneous and isotropic 
Friedmann-Robertson- Walker (FRW) space-time. How- 
ever, in this theory there are subtleties arising due to 
the freedom of gauge, i.e. choosing the correspondence 
between perturbed and background quantities 0, 0] . 

Another approach is Newtonian cosmology, see e.g. [I]- 
One has to keep in mind that Newtonian theory is wrong 
in that it assumes instantaneous gravitational interaction 
and infinite speed of light. However, because of its sim- 
plicity compared to the relativistic theory, the Newtonian 
equations are the preferred choice in cosmological N-body 
simulations, which are used to test our understanding of 
structure formation. 

But the question is, how reliable are these simulations, 
since they use Newtonian rather than relativistic equa- 
tions. Here we ask, how reliable is Newtonian cosmology 
on large scales? 
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In order to answer this question, people have so far 
followed mainly two different strategies. One strategy 
is to choose a gauge and then calculate the relativistic 
corrections appearing in this gauge. It turns out that in 
any gauge that is not the synchronous gauge there are 
relativistic corrections appearing to the density contrast 
on large scales, see e.g. |9l. [Toj. 

The other strategy is based on the use of gauge- 
invariant variables. Gauge-invariant quantities were first 
introduced by Bardeen in Q (see also Q for a more re- 
cent review) and are the basis for the work of Hwang 
and Noh, who found an exact correspondence up to sec- 
ond order between Newtonian and relativistic cosmology 

The relativistic-Newtonian correspondence up to sec- 
ond order found in |11| holds in the following way: The 
Newtonian density contrast is identified with a relativis- 
tic gauge-invariant expression that reduces to the density 
contrast on comoving hypersurfaces. However, the New- 
tonian gravitational and velocity-potentials are identified 
with gauge-invariant expressions which reduce to the rel- 
ativistic potentials on zero-shear hypersurfaces. Using 
the same strategy, Hwang and Noh also found a corre- 
spondence between relativistic cosmological perturbation 
theory and the first order post-Newtonian approximation 
(IPN) [13. Again, this correspondence holds if one iden- 
tifies the IPN perturbation fields with relativistic gauge- 
invariant quantities that are interpreted on different hy- 
persurfaces. 

However, one needs to keep in mind that any phys- 
ical problem is described by a set of equations of mo- 
tion and appropriate initial conditions. The initial condi- 
tions must be specified on a spatial Cauchy hypersurface, 
which in the context of cosmological perturbation theory 
corresponds to a particular foliation of space-time, i.e. to 
a hypothetical observer who is able to determine physi- 
cal quantities on a spatial hypersurface. The relativistic- 
Newtonian correspondence mixes the quantities defined 
on different spatial hypersurfaces and thus no hypothet- 
ical observer in the Einsteinian world could actually de- 



2 



tcrminc these combined quantities. 

Another recent approach from Green and Wald uses a 
whole new framework to investigate the problem [S] ■ 
Using their framework, the authors construct a dictio- 
nary between Newtonian solutions and relativistic solu- 
tions. However, in order to reconstruct the exact full 
relativistic structure out of the Newtonian solution, ad- 
ditional differential equations need to be solved (e.g. eq. 
(2.44) in ref. [lEl)- A similar dictionary is proposed by 
Chisari and Zaldarriaga [l^ . They show how relativistic 
corrections can be absorbed into initial data of simula- 
tions and how to modify the Newtonian coordinates to be 
able to compare them to relativistic coordinates. Both 
strategies seem to result in extra input and extra numer- 
ical work on top of the Newtonian calculation. 

In this work we give a quantitative estimate of when 
the relativistic corrections to the Newtonian solutions 
are important. In order to do so, we first construct a 
gauge in which both the density contrast and the veloc- 
ity perturbations coincide with Newtonian theory. Then, 
the significance of the relativistic corrections to Newto- 
nian theory can be evaluated by considering the local 
modification of the Hubble parameter, which appears in 
this gauge, and comparing it to second order Newtonian 
quantities. 

Let us note that the question of this work is different 
from the question what a real observer can see in the 
Universe. The latter has been recently discussed in (TtI - 

Throughout this work we restrict our attention to an 
Einstein-de Sitter background cosmology, because it pro- 
vides an accurate description for most of the history of 
the Universe and has simple expressions for the scale 
factor (a oc t^, where r denotes conformal time) and 
the conformal Hubble parameter {H — 2/t). We adopt 
a Hubble constant of Hq = 71 km/s/Mpc, in agree- 
ment with measurements of the Wilkinson Microwave 
Anisotropy Probe (WMAP) [2l|. 

The outline is as follows: In section II we review New- 
tonian cosmological perturbation theory and present the 
equations up to second order. In section III we describe 
relativistic cosmological perturbation theory and give an 
overview over solutions in the six most common gauges. 
In particular, we show how the shapes of power spectra 
vary from gauge to gauge. We find that none of these 
gauges gives both matter- and velocity spectra in cor- 
respondence with Newtonian theory. In section IV we 
introduce a new gauge with the property that it coin- 
cides with Newtonian theory on all scales regarding the 
density contrast and the velocity perturbation, the New- 
tonian matter gauge. Using the strategy described above, 
we conclude that relativistic corrections are more impor- 
tant than second order Newtonian effects on scales larger 
than 10 Mpc. 



II. NEWTONIAN COSMOLOGICAL 
PERTURBATION THEORY 

First order 

Newtonian cosmological perturbation theory results 
in perturbed Newtonian equations (continuity equation, 
Euler equation and Poisson equation) on an expanding 
FRW background. In this model the linear equations can 
be written in Fourier space Q: 



5W' + ik-v(i) =0, 



(1) ^ in'6^'\ 



(1) 

(2) 
(3) 



where S^^"^ is the first order density contrast, v'^-'^^ is the 
first order velocity perturbation, 0'^^ is the first order 
gravitational potential and k is the comoving wavenum- 
ber. The linear equations can be decoupled into scalar 
and vector parts. It turns out that in an Einstein-de Sit- 
ter model the vector contribution to v^^^ is decaying and 
hence can be neglected (see e.g. HH). 

The remaining growing solutions for the perturbations 
are: 



.(1) 



2ik 



0W. 



(4) 
(5) 
(6) 



We assume a Harrison-Zel'dovich spectrum for the pri- 
mordial curvature perturbation field, with an amplitude 
A = 4.88 • 10~^, which corresponds to measurements of 
WMAP HH. It follows that the initial condition for ^(^^ 
can be written as (see e.g. [2^): 



|^(i)| = -^(27r2)i/2fc-3/2T(fc), 



(7) 



where T{k) is the transfer function, which accounts for 
suppressed growth of structure during the radiation dom- 
inated epoch. We adopt the form of T{k) from [23] for a 
pure cold dark matter model (i.e., without any baryonic 
matter) . 



Second order 

At second order the Newtonian equations are more 
complicated. For a detailed derivation, see [l^. In 
Fourier space they read 



v(2)'+7{v(2)+ik</)(2) 



ik 



(27r)= 

ik 
(2^ 



■ ((5(i'*v(i)), (8) 



(v(i)*v(i)), (9) 



(10) 
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where a star denotes a convolution. Note that the struc- 
ture of the second order equations is very similar to that 
of the first order equations. However, there are additional 
source terms, consisting of convolutions of first order per- 
turbations, which we have written on the r.h.s. of each 
equation. Due to these source terms, the second order 
perturbations grow faster than the first order perturba- 
tions. Solutions can be obtained numerically. Plots of 
the first and second order density contrast are shown be- 
low in figure [5J where we compare them to a relativistic 
quantity, the local variation of the Hubble parameter. 



III. RELATIVISTIC COSMOLOGICAL 
PERTURBATION THEORY 

In relativistic cosmological perturbation theory we 
consider small perturbations of the homogeneous and 
isotropic metric and the energy-momentum tensor of a 
perfect fluid. The evolution of these perturbations follows 

I 



from the covariant conservation of the energy-momentum 
tensor and Einstein's field equations. We define the per- 
turbed metric 



-1 — 20 W;i 

w.,, (1 - 2V')5y -I- 2h.. 



(11) 



We restrict our attention to scalar perturbations at linear 
order. There are hence 4 degrees of freedom in the metric 
perturbation, corresponding to 0, ip, w and h. However, 
two of them can be fixed due to the gauge freedom. This 
corresponds to choosing and ^ in the gauge transfor- 
mation 



x + 



(12) 



where t is the conformal time and x are comoving coor- 
dinates. In the following we will discuss the advantages 
of choosing w = as a gauge condition. The proper time 
between two events along a worldline is given by 



2<i)-2w~ - [(1 

CtT 



2i))5,j+2h,,j\ — — . 



(13) 
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Now we expand the integrand and keep only terms up to 
0{gvi), where g S {(j>, ip, w, h, Vi}, which gives 



— V • Vw 

2 



(14) 



The proper time is extremal if and only if the Lagrangian 



L = a{l- 
equation 



satisfies the Euler- Lagrange 



d dL _dL 
dr (9v 9x 

This gives in first order 

1 rf[a(v + Vw)] _ 
a dr 



-V0. 



(15) 



(16) 



It is convenient to choose w = 0, so that the Euler- 
Lagrange equation coincides exactly with Newton's law 
for the motion of particles in an expanding universe. 



1 d(av) 
a dr 



(17) 



Thus, in a gauge with w = we can identify the metric 
perturbation 4> with the Newtonian gravitational poten- 
tial. This is the reason why it is common to identify the 
Bardeen potential <&, which is the same as (j) in the lon- 
gitudinal gauge (where w = ft- = 0), as the gravitational 
potential (the exact definition of $ is given below). In 
physical terms, w = implies that comoving observers 



do not change their coordinates between constant-time 
hypersurfaces. 

After choosing w = there is still another degree of 
freedom, which can be fixed by a second gauge condition, 
e.g. 

• Uniform curvature gauge (UC). In this gauge we 
set w = i/) = 0. It follows that the intrinsic Ricci 
curvature, which is given by '"^^i? = ^A^/;, vanishes 
in this gauge 0- 

• Longitudinal gauge (L). In this gauge we set w = 
/i = 0. It is equivalent to the zero shear gauge, 
defined by if = x = 0, where 



X = a{h' — w) 



(18) 



generates the traceless part of the extrinsic curva- 
ture tensor and can be interpreted as shear in the 
normal worldlines Another common name for 
this gauge is conformal Newtonian gauge. How- 
ever, we will not use this name since it is not the 
only gauge where the particles move according to 
Newton's law of motion, as we have shown above. 
Bertschingcr's Poisson gauge reduces to the longi- 
tudinal gauge in the scalar sector @. 

The longitudinal gauge is equivalent to the use of 
the Bardeen potentials 0, 



$ = + i [(w - /i')a]' . 
a 

^ = ^P-niw-h'). 



(19) 
(20) 
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FIG. 1. Matter power spectrum [Pss = matter-velocity power spectrum {PvS = |v(5|) and velocity power spectrum 

(Pvv = |vp) according to relativistic perturbation theory in different gauges (defined in the text) and according to Newtonian 
theory (labelled by "N"), for 2; = (left panels) and z = 100 (right panels). 
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TABLE I. Overview - growing modes of first order scalar perturbations for dust models according to relativistic cosmological 
perturbation theory in longitudinal (L), uniform curvature (UC), spatially Euclidean (SE), synchronous (S), comoving (C), 
uniform density (UD), and uniform expansion (UE) gauge. The quantities are arranged in the following way: The upper group 
contains perturbations in the metric tensor, the middle group contains perturbations in the energy-momentum tensor (peculiar 
velocity potential and density contrast) and the lower group contains derived geometrical quantities (intrinsic curvature, shear, 
and expansion rate). All solutions are expressed as function of Bardeen's metric potential $. 



• Uniform expansion gauge (UE). In this gauge we 
set w ~ K ^ 0, where 

ac=-(V'' + H0)-4ax (21) 

is the perturbation of the trace of the extrinsic cur- 
vature tensor 0]. Since the full trace is given by 
K = —3H + K, where H is the Hubble expansion 
rate, we can identify SH = — k/3 as the perturba- 
tion in the expansion rate. 

• Uniform density gauge (UD). In this gauge we set 
w = d ^ 0. There are hence no density perturba- 
tions. 

• Synchronous gauge (S). In this gauge we set w = 
(j) = 0. It follows that observers at different places 
have synchronous clocks, since affects the time- 
time component of the metric tensor. There is a 
residual gauge freedom, which can easily be fixed 
by specifying that the peculiar velocity vanishes at 
a particular moment of time (see e.g. Q.) For a 
matter dominated universe this coincides with the 
comoving gauge (see below). 

• Comoving gauge (C). In this gauge we set w = v = 
0, where v is the velocity potential, so that v = Vf . 
There are hence no peculiar velocities in this gauge. 
It turns out that during matter domination the co- 
moving gauge and the synchronous gauge coincide 
[2^ . Hence, we will treat them as one gauge in the 
following. 

• An interesting gauge with w ^ is the spatially 
Euclidean gauge (SE). Here we set ft, = -0 = 0, so 
that the spatial part of the metric has Euclidean 
geometry. 



The relativistic equations of motion for the fluctua- 
tions follow from the perturbed energy-momentum con- 
servation equation, Tj(^^ = 0, and the Einstein equation, 
Gfii, = SttGT^^, where Gf^, denotes the Einstein tensor 
and T^i, is the energy-momentum tensor. In longitudinal 
gauge, covariant conservation of energy and momentum 
gives for dust (see e.g. Q) 

(5^i^L + AwGi.L - 3*' = 0, (22) 
Vl'GI.L + ^Vl-GLL = -V$, (23) 



and the Einstein equation gives 





3-H^$ + 3^^-' - 


-AvI/ = -^H25gi.l,(24) 






= -^H2i;gi,l,(25) 


3*" + 


3H($' + 2^'') + A(* 


- $) = 0, (26) 




(VI/™ 


= 0, (27) 


where we 


have introduced the 


gauge-invariant density 



contrast, 



^Gi,L ^S- 3n{w - h'), (28) 

as well as the gauge-invariant velocity potential, 

WGi.L = V + h'. (29) 

Thus, the equations are written in gauge-invariant form. 
Note that the fields $, (5gi,l, wgi.l reduce simply to 
0L, (^L, wl, respectively, in the longitudinal gauge. This 
tells us how to interpret this set of gauge-invariant fields: 
They measure the perturbations on zero-shear hypersur- 
faces. 
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In longitudinal gauge, the growing solutions are: 



IV. NEWTONIAN MATTER GAUGE 



V'L = 0L = "J* = constant, 
2 $ 



Vh 



(30) 
(31) 

(32) 



The solutions in other gauges can be obtained using 
the following gauge transformations 0: 



X - X - 

2 a 



(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
(39) 

(40) 



The solutions in different gauges are given in table H 

Figure [1] shows plots of the the matter power spec- 
trum, Pss = I^P, the matter- velocity power spectrum, 
PvS = and the velocity power spectrum, Pw = |vp. 
Each spectrum is shown for different gauges and at two 
different redshifts: z — 100, because this is where typi- 
cal numerical simulations begin, and z = 0, because this 
is where typical numerical simulations end. The initial 
condition for $ is given by equation ([7]). This is possi- 
ble because we can identify $ = 0l as the Newtonian 
gravitational potential, as seen from (|17p . 

Figure [T] demonstrates that the relativistic matter 
power spectrum coincides with the Newtonian one on 
all scales in synclironous/comoving gauge, however there 
are no relativistic velocity perturbations. For the longi- 
tudinal gauge and the spatially Euclidean gauge the rela- 
tivistic velocity power spectrum coincides with the New- 
tonian one on all scales, but the matter power spectra do 
not match on superhorizon scales. The other gauges give 
non-Newtonian results in both spectra. None of these 
gauges gives a correspondence to Newtonian theory both 
in the matter- and the velocity power spectrum on all 
scales. 

Note that the claims about the Newtonian-relativistic 
correspondence found in literature are based on either the 
use of appropriate dictionaries [isl [l6| or the appropriate 
mixing of gauge- invariant variables jl2|. We stress that 
the findings in this section are solutions as measured by 
one hypothetical observer, i.e. in one gauge. It seems to 
us that, in order to estimate the relativistic corrections to 
Newtonian cosmology, first of all we must seek a gauge in 
which both the density contrast and the matter velocity 
coincide. 



By choosing and ^ adequately, it is possible to con- 
struct a gauge in which both the density contrast and the 
peculiar velocity coincide with their Newtonian counter- 
parts on all scales. We call this gauge the Newtonian 
matter gauge (NM), as all variables associated with the 
state of the matter (dust) agree with their Newtonian val- 
ues. In order to construct it, we start in the longitudinal 
gauge, where the peculiar velocity already coincides with 
the Newtonian one. We then choose = |^ so that 



3NM 



(41) 



Note that the (0, 0)-component of the metric perturba- 
tion vanishes in this gauge. 



0' 



0, 



(42) 



which is not a problem because the (0, i)-part of the met- 
ric perturbation is not zero in this gauge, 

2<j> 

WNM=wi^+e-C'=C' = ^, (43) 

which means that we cannot identify 0nm as the Newto- 
nian gravitational potential, as we discussed earlier. 

It is illuminating to evaluate (|16p in the Newtonian 
matter gauge. Since 0nm = 0, we have 



Id 

- — [a(vNM + Vwnm)J = 0. 
a uT 



Using 



Id, . ld/2$\ 1^3* 



(44) 



--(au;NM) = -^(^a-j=-,-r^- = a> (45) 

we recover Newton's equation of motion. This is not a 
surprise because we have chosen the gauge to give matter 
velocities in agreement with Newtonian theory. 

In the Newtonian matter gauge we find a non- vanishing 
intrinsic curvature. 



3 a-' 



(46) 



as well as a non-vanishing expansion rate perturbation 

2 fc2 



and non- vanishing shear, 

Xnm = 



3 an 



2 a 




(47) 



(48) 



All these quantities do not appear in Newtonian cosmol- 
ogy. In particular, we define the local modification of the 
Hubble expansion rate. 



OH 



(SH) 



NM 



H 



3H ' 



(49) 
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FIG. 2. Density contrast in first and second order Newtonian cosmological perturbation theory, compared to the quantity 5h 
(defined in the text), which represents the local modification of the Hubble rate and is a purely relativistic effect. The prefactor 
is chosen to make the quantities dimensionless. The different panels show the behaviour on different scales: on comoving length 
scales larger than 10 Mpc the relativistic modification is always larger than second order Newtonian effects, while on scales 
below 1 Mpc the relativistic modifications are small compared to second order Newtonian effects. These results hold for a 
spatially flat dust universe. 



Since the density contrast and the peculiar velocities 
coincide with Newtonian theory in this gauge, we can 
use 5h to estimate the magnitude of relativistic modifi- 
cations to the Newtonian theory. In Figure [2] we show 
the density contrast at first and second order Newtonian 
cosmological perturbation theory, as well as 5h on differ- 
ent scales. We see that on comoving length scales smaller 
than 1 Mpc the second order Newtonian effects are al- 
ways larger than the relativistic modifications. However, 
on comoving length scales larger than 10 Mpc the rela- 
tivistic modification is more important than the second 
Newtonian order. This can have a drastic effect on the 
reliability of Newtonian N-body simulations of large cos- 
mological volumes. Thus, Newtonian and general rela- 
tivistic models of cosmology could coincide in the matter 



and velocity power spectra for all times, but would then 
disagree in their Hubble diagrams and thus in their con- 
clusions on the expansion history of the Universe. 

In order to make a quantitative statement about the 
validity of Newtonian theory we define two characteristic 
redshifts. Let znl be the redshift at which the perturba- 
tion on a given scale goes nonlinear, which is, according 
to the spherical collapse model [l^, given by 



1, 



(50) 



where V 



(1) _ 



r|i5^"'^^P is the dimensionless matter 



5 — 271^ 

power spectrum. Perturbation theory can only be ap- 
pUed up to this redshift. For comparison, let znn be the 
redshift at which the relativistic modifications become 
as important as the second order Newtonian corrections. 
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FIG. 3. The local modification of the Hubble parameter, 5h, 
plotted against the comoving wavenumber k. From top to 
bottom, the lines denote growing redshifts of 0, 2, 5, 10, 20 
and 50, respectively. Solid lines mark the regime of validity of 
linear perturbation theory. Our results do not strikly apply 
in the nonlinear regime (dotted lines). 
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TABLE II. 2NL and znn (defined in the text) for different co- 
moving wavenumbers. For k > 0.6 Mpc^^ we have znl > znn, 
which means that these scales go nonlinear before relativistic 
effects become more important than second-order Newtonian 
effects. However, for k < O.BMpc"^ we find znl < ^nn, so 
that relativistic corrections become more important than the 
second order Newtonian terms before these scales go nonlin- 
ear. 

that is, 

(5^^^(fc,ZNN) '-^ <5_fi-(fc, 2nn)- (51) 

In table HI] we show an overview of these characteristic 
redshifts for different comoving scales, which can be read 
off from figureHl using the relation a = l/(l-|-z) between 
scale factor and redshift. A discussion follows in the next 
section. 

In figure [3] we plot the dimcnsionless band power of 
5h against the comoving wavenumber k for different red- 
shifts. It can be seen that the relativistic corrections to 
the Hubble rate on the cluster scale (fc ^ 0.1 Mpc^^) is 
of the order of 1% at a redshift of z ~ 2 and 10% today. 

Note that relativistic corrections on the Gpc-scale are 
always below 1%. However, this result should not be in- 
terpreted too fast as a confirmation of Newtonian theory 
on these very large scales. In fact, the density contrast 
itself is very small on the these scales. In particular, 



the fraction Sh/S^^'' = a^/3 is independent of scale, so 
that the relativistic corrections become comparable to 
the density contrast today on all scales that still behave 
linear. 



V. DISCUSSION 

The theory of structure formation gives the transition 
between the early homogeneous Universe and the Uni- 
verse we observe today. In this work we have studied 
two different theories that describe this transition: rela- 
tivistic and Newtonian cosmological perturbation theory. 
Although the former is the physically more correct tool, 
the latter is, due to its simplicity, the preferred choice in 
cosmological N-body simulations, with which we test our 
understanding of structure formation. Fixing the back- 
ground universe to be an Einstein-de Sitter model, we 
have calculated the first and second order perturbations 
in Newtonian cosmological perturbation theory and the 
first order perturbations in relativistic cosmological per- 
turbation theory. 

We have shown how the power spectra for mat- 
ter and velocity perturbations behave according to rel- 
ativistic cosmological perturbation theory in different 
gauges and according to Newtonian theory. In the syn- 
chronous/comoving gauge the density perturbations be- 
have exactly like in Newtonian theory, however there are 
no velocity perturbations. In the longitudinal gauge the 
velocity perturbations behave like in Newtonian theory, 
but the density contrast differs significantly on superhori- 
zon scales. 

In the newly defined Newtonian matter gauge there 
are no relativistic corrections at all to the Newtonian 
matter- and velocity power spectra. However, there are 
other quantities present which do not appear in Newto- 
nian cosmology: shear, intrinsic curvature and pertur- 
bations in the expansion rate. In particular, we have 
introduced the local modification of the expansion rate, 
Sh , and compared it to second order terms in Newtonian 
theory in order to evaluate at which scales Newtonian 
theory represents a good approximation. 

Clearly, Newtonian theory is a good approximation 
when the scale goes nonlinear before relativistic modi- 
fications can become more important than second order 
Newtonian corrections, i.e. for znl > ^nn- This is true 
for k > 0.6Mpc~^, see table HIl For znl < ^^nn, relativis- 
tic modifications become more important than second or- 
der Newtonian corrections before the corresponding scale 
goes nonlinear. This is the case for k < 0.6 Mpc^^. 

Here, we only focussed on one relativistic effect, the lo- 
cal variation of the Hubble parameter. There are however 
also other relativistic quantities that Newtonian simula- 
tions do not "see", namely intrinsic curvature and shear. 
When the fluctuations in the Hubble expansion rate be- 
come relevant, the tracing of light rays through a Newto- 
nian simulation also must start to deviate from what the 
full theory would predict. This has to be kept in mind 
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when using huge volume simulations in quantitative com- 
parison with real data. 

We stress that our results do not contradict the claims 
about the exact Newtonian-relativistic correspondence of 
equations of motions written in terms of Bardeen vari- 
ables in different foliations found in literature, e.g. by 
Hwang and Noh [ll|, Ell. We have calculated the rel- 
ativistic corrections as they appear in one specific and 
fixed foliation of space-time, i.e. as measured by one 
specific hypothetical observer in the linearly perturbed 
world, while the work by Hwang and Noh is based on 
the combination of gauge-invariant variables that can be 
interpreted on different hypersurfaces. 

We are now in a position to be able to quantify the 
reliability of Newtonian simulations at large scales. We 
can interpret existing Newtonian simulations (i.e. simu- 
lations that have not been corrected or modified along the 
lines proposed in [l5l.[l6j) to correspond to simulations in 
the Newtonian matter gauge. Then the simulated density 
contrast and peculiar velocities are correct (by definition) 
and other observables like the Hubble expansion rate are 
modeled very well on scales below 10 Mpc. However, on 
larger scales relativistic corrections are more important 
than Newtonian non-linear effects. 

Thus we conclude, reliable predictions on supercluster 
scales, void scales, the baryon acoustic oscillation scale 
and beyond can only be based on general relativistic equa- 
tions. Newtonian simulations are good to provide us with 
a qualitative picture at scales above 10 Mpc, but a mea- 
surement of cosmological parameters at better than 10 
per cent level cannot rely on them. Our findings do not 



question that Newtonian simulations remain to be im- 
portant and extremely valuable at the cluster scale and 
below, i.e. in the deeply non-linear regime. 

In Newtonian cosmology, the so-called redshift space 
distortions [l^ are due to the peculiar motion of galax- 
ies, whereas in the relativistic formulation in Newtonian 
matter gauge, they would receive contributions from the 
peculiar motion and from space-time metric perturba- 
tions, most importantly the local variation of the cosmic 
expansion rate Sh- This clearly demonstrates that the 
differences between Newtonian and relativistic dust cos- 
mologies should be observable in high fidelity galaxy red- 
shift surveys. A quantitative analysis of these differences 
requires further investigations. 

Let us finally remark that also exact theorems support 
our point of view that general relativistic and Newtonian 
dust cosmology are inequivalent. It has been shown by 
Ellis that there are no shear- free dust solutions of the Ein- 
stein field equations that both expand and rotate (" Dust 
Shear- Free Theorem", see [1^). However, solutions of 
this kind do exist in Newtonian cosmology, as was shown 
by Narlikar This alone demonstrates that the two 

theories are inequivalent. 
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